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Abstract. By means of the time-dependent density matrix renormalization group 
algorithm we study the zero-temperature dynamics of the Von Neumann entropy of a 
block of spins in a Heisenberg chain after a sudden quench in the anisotropy parameter. 
In the absence of any disorder the block entropy increases linearly with time and then 
saturates. We analyze the velocity of propagation of the entanglement as a function of 
the initial and final anisotropies and compare, wherever possible, our results with those 
obtained by means of Conformal Field Theory. In the disordered case we find a slower 
(logarithmic) evolution which may signals the onset of entanglement localization. 
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1. Introduction 

The recent interest in aspects common to quantum information and condensed matter [1] 
has prompted a florishing activity at the border of these disciplines, quite distinct 
until few years ago. It is hard to list all the numerous problems addressed so far. 
Here we only mention the study of entanglement in quantum critical systems (see for 
example [2-23] and references therein), and the recent advances relating to the density 
matrix renormalization group (DMRG) [24-26], as they are tightly connected to the 
present work. In fact in this paper we apply the recently developed time dependent 
DMRG (t-DMRG) [24-26] to the study of the dynamics of entanglement entropy in 
Heisenberg chains. 

Among all the various ways to quantify entanglement here we consider the block 
entropy which has recently analyzed in several different situations (see e.g. [4,9,11,17- 
19,27-30]). It has been demonstrated [4,9,31] that the entropy of a block of spins St 
in the ground state of a spin chains is very sensitive to its critical properties. In the 
case of a block with one boundary with the rest of the chain Si, for £ 3> 1, diverges at 
the critical point as Si = (c/6) log 2 £ where c is the central charge of the corresponding 
conformal field theory (CFT) of the model considered. In contrast in non critical systems 
the entropy saturates to the finite value St = (c/6) log 2 £, with £ 1 the correlation 
length [9]. In the presence of quenched disorder the properties of the block entropy for 
critical chains remain remarkably universal. For the models analyzed in Ref. [32] using 
real-space renormalization group, it still diverges logarithmically with a "renormalized" 
central charge c e g- = cln2 [32], where c is the central charge of the corresponding 
pure model. It has been conjectured that such an effective central charge characterizes 
generically the critical systems with quenched disorder [32,33]. 

The interest in the properties of entanglement in condensed matter has also 
extended to understand its dynamical behaviour. Like for the case of propagation of 
excitations in condensed media, it recently became of interest to know how entanglement 
could propagate through spin chains. This question was studied by looking at two- 
particle [34-37] and many-particle entanglement [38,39]. The dynamics of entanglement 
was studied either by preparing the system in a state (not an eigenstate) with all the 
entanglement localized in a given part of the chain or after a sudden quench of some of 
the couplings of the model Hamiltonian. 

In this paper we consider static and dynamical properties of the entropy of a block 
of spins for an anisotropic Heisenberg chain both in the clean and in the disordered 
case. We investigate the problem by means of DMRG [40] for the static part [41] and 
its time-dependent version [25, 26] for the dynamical evolution. The static entropy 
of the Heisenberg model was considered numerically in Ref. [4, 19, 29, 42, 43] and some 
exact results are known in the isotropic ferromagnetic limit [44] and in the XX and Ising 
limits [6,13,18,19]. In this work we consider again this case for completeness. In addition 
it serves as an important check for the other cases considered (see Sections I3.1| 13 .2\ 14. 1| 
I4.2JI for which these are, to the best of our knowledge, the first numerical calculations 
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available for the Heisenberg model. The results of the dynamics of an ordered chains 
are compared with the CFT results of Ref. [38]. In the static disordered case we confirm 
the prediction of [32]. Section l4~2l is devoted to the dynamics of a disordered Heisenberg 
chain. Here we predict a slow evolution for the dynamics of entanglement which hints 
at some sort of entanglement localization. 



2. The model 



The model that we consider is a spin chain of length iV with open boundary conditions 
described by the Heisenberg Hamiltonian [45] 

N-l 
i=l 

where c l x y z are the Pauli matrix operators relative to the ith spin; Jj are the coupling 
constants that we assume time independent but possibly space-dependent in the random 
case; finally A is the anisotropy parameter. The homogeneous chain is critical for 
— 1 < A < 1. The corresponding Conformal Field Theory (CFT) is described using a 
central charge c = 1. 

We study the properties of the Von Neumann entropy Si of a block containing the 
first £ spins. The block is described by the reduced density matrix p t = Tr^p and Se 
is defined as: 

St = -Tr £ (pt log 2 p t ) (2) 

We used the t-DMRG algorithm with a second order Trotter expansion of H as described 
in Refs. [25,26]. We checked the precision of the numerics by comparing the results with 
the case A = where Eq. (jlj is mapped onto a free fermion model [45]. 

We first present the result for the homogeneous case and then that for the disordered 

one. 



3. Homogeneous chain 

3.1. Ground state properties 
The CFT prediction for a block of length 



i log2 
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in an open chain of total length N is [9] 
+ A (3) 



where A is a non-universal constant related to the analogous one in the system with 
periodic boundary condition [9, 29]. As already mentioned the static properties of 
the Heisenberg model was considered numerically in Ref. [4, 19, 29, 42] and here for 
completeness we present the result of our DMRG calculations. They are shown in Fig. 
^ The main plot shows the difference in the entropy for a critical and a non critical value 
of A and the comparison to the CFT prediction. For large I (i.e. i > 10) where Eq. © 
is expected to work, the agreement is very good. In the lower inset we show our results 
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for the central charge, obtained through a fit of numerical data with Eq. (jSJ), for different 
values of A in the critical interval. The fitted c shows small non-universal variations 
with A which decrease increasing N and are expected to vanish in the thermodynamic 
limit. This is shown in the top inset for the worst case A = 0.5. The extrapolated value 
of the central charge (upper inset) has been obtained fitting c(l/N) with a quadratic 
polynomial and taking the limit 1/N — > 0. The result is c = 1.01 ± 0.05. Note that 




i i i i i ' ■ ' i i i i i i i i i 

1 10 £ 100 

Figure 1. The Block entropy Si for N — 200 for a critical value A = 0.0 (circles) and 
non-critical value A = 1.8 (squares) and m = 120. The critical data compared with 
the CFT prediction Eq. © (dashed line). Lower inset: central charge extrapolated by 
fitting the numerical data Si with Eq. for different values of A. The data are for 
N = 1000 and m = 120. Upper inset: scaling of c extrapolated as a function of 1/N 
for the worst case A = 0.5 and compared to a quadratic fit (dashed line). 

since we are using open boundary conditions the entropy oscillates with the parity of the 
block. This can be explained with a simple argument. Because of the XXZ interaction, 
neighbor spins tend to form spin singlets, so the oscillations in the entropy of block 
reflect the breaking or not of one of these singlets. Such an alternating behaviour has 
been observed also in Ref. [43]. Though this behavior disappears increasing N we fitted 
numerical data for even I because the convergence is faster. 

3.2. Dynamical behaviour 

To date the only results obtained for the evolution of block entropy after a quench 
have been obtained in Ref. [38]. By means of CFT it was shown that a quench of the 
system from a non critical to a critical point leads the block entropy to increase in 
time until a saturation point is obtained. For periodic boundary conditions, the time at 
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which the entropy saturates is given by t* = £/(2v) where v is the spin wave velocity: 
v = dE^I dk\k = Q. This phenomenon has a simple interpretation in terms of quasiparticles 
excitations emitted from the initial state at t = and freely propagating with velocity 
v. This explanation holds even for non critical systems, as it has been confirmed by 
the exact solution of the Ising chain dynamics [38]. However, since in lattice models 
there are particles moving slower than v, after t* the entropy does not saturate abruptly, 
but is a slowly increasing function of the time. We now consider the dynamics of the 
Heisenberg model with open boundary conditions, and we will find that t* = £/v. At 
the end of this section we will interpret this result in terms of quasiparticles and we will 
derive it within CFT. 



3 
2.5 
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Figure 2. Evolution of the entropy Sq with various quenches. Ao = 1.5 while 
Ai = 0.0,0.2,0.4,0.6,0.8 as a function of u(Ai)i. Inset: initial slope value of Sq 
as a function of Ai and comparison to a linear fit with slope —0.85 ± 0.02 (dashed 
line). 




In the case of the Heisenberg Hamiltonian in the critical regime the spin wave phase 
velocity is given by v(A) = 2Jtt sin6/6 with cos 6* = A [46]. The initial state of the 
system is chosen as the ground state of Hamiltonian Eq. with A = A > 1, after 
a quench the system evolves with same Hamiltonian but with a different anisotropy 
A = Ai G [0; 1]. In the simulations we considered chains with N = 50, a Trotter slicing 
JSt = 5 ■ and a truncated Hilbert space of m — 200. The block was chosen to be 
of 6 sites which is large enough, as we show, to confirm the CFT prediction. We have 
checked convergence with m and St. For the special case A = +oo and Ai = we 
compare our data with the exact result obtained diagonalizing the XX model as shown 
in Appendix A 
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Figure 3. Evolution of the entropy Sq with various quenches. Ao = 1.2, 1.5, 3.0, oo 
while Ai = 0.0 as a function of v(Ai)t and shifted so to coincide in t = 0. For A = oo 
we show also the exact result obtained by diagonalization (circles ). Inset: initial slope 
value of as a function of Ai. 

Various quenches have been considered. First of all we present the results with A 
fixed and Ai variable shown in Fig. El As in Ref. [38] the entropy grows linearly in 
time S ~ sJt and finally saturates. At t = t* = £/v, the entropy does not saturate 
abruptly. As discussed for the Ising model, this is due to slow quasiparticles [38]. The 
inset of Fig. El shows the initial slope as a function of Ai. The dashed line is a linear fit 
with slope 0.85 on the data except A x = 0. In Fig. El the time evolution of the entropy 
is shown for quenches starting from different values of Ao and ending with the same 
Ai = 0. The case in which both A and Aj are varied (keeping their difference fixed, 
Ao — Ai = 1.5) is shown in Fig. |U 

The behaviour of the slope s (the entropy for t < t* goes as S ~ sJt ) for I = 20 
as a function of A and A x can be fitted for a large portion of values by the law 

s = 1.50A - 0.84A! - 0.90 . 

In Fig. El the fit and the numerical data are compared. Following the simple model 
introduced in [38], the initial slope of the entropy increase is a non-trivial function of 
both A and A x as it depends both on the crossection for producing the quasiparticles 
and their velocity. 

According to the argument reported in Ref. [38], the block entropy is proportional 
to the numbers of pairs of entangled quasiparticles, emitted from any point at the time 
t — 0, of which one reaches the block and the other reaches the rest of the system. If 



Entanglement Entropy dynamics in Heisenberg chains 



7 





Figure 5. Derivative s of the entropy with respect to time as a function of Ao and 
Ai (surface) and compared to a linear fit (grid lines): s — 1.50Ao — 0.84Ai — 0.90. 
The deviation with the linear fit is less than 10% for (A , Ai) e [1.4,2.4] x [0.4,0.8]. 



these particles are all moving at speed v, the entropy of a block in a chain infinite in 
both directions is at first linear in t until t = t* = l/(2v) and then saturates to a value 
proportional to I. Clearly the presence of boundary condition at x — 0, changes this 
result, since it acts as a wall and particles reaching it cannot propagate any longer in a 
straight line. If we assume that the wall at x = is reflecting, it works like a mirror for 
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the quasiparticles. The resulting "effective length" the block is 2£ and the saturation 
time t* = 2£/(2v), in agreement with what we have found numerically. On the other 
hand a completely adsorbing wall does not change the "effective length" the block and 
the resulting saturation time is the same as for periodic conditions. 

It is very simple to generalize the CFT calculation of Ref. [38] to the case of a 
boundary condition that is of the same kind of the initial condition (e.g. fixing the 
spins in the same direction both at initial time and at x = 0, or having a disordered 
initial state and free boundary condition at x = 0). In this shown in [9], Trp^ 

transforms like a one-point function. The imaginary-time half-strip of width 2e (in the 
notation of Ref. [38] to which we refer the reader for details) is mapped on the half- 
plane by the conformal transformation z(w) = sin(nw / (2e)) (more details about this 
calculation can be found in Ref. [47] where it has been performed in a different context). 
Continuing the one-point function obtained to real time and performing the "replica 
trick" [9], we easily obtain t* = £, in agreement with the quasiparticles interpretation 
and our numerical results. 

4. Disordered chain 

4-1. Ground state properties 

The striking prediction of Refael and Moore [32] that the block entropy of random 
chains has a logarithmic divergence with an effective (universal) central charge has been 
confirmed by numerical calculations on the XX model [33]. Here we extend these 
results to the model defined in Eq. (JTJ) where the randomness is introduced through the 
couplings Jj chosen to be random numbers with a uniform distribution in the interval 
[0; 1] J. According to Ref. [32], the block entropy should scale with an effective central 
charge c e g = In 2 ~ 0.69. In Fig. El we show for N = 50 the entropy of a block of 
length I for the random model for A = 1 averaged over 10 4 different configurations of 
disorder. From the scaling of the entropy as a function of i we can extract the effective 
central charge for this model. The result is cgt = 0.67 ± 0.05 in very good agreement 
with Refael and Moore proposal. 

4-2. Dynamical behaviour 

No analytic results are known so far for the dynamics of block entropy in the case 
of disordered systems. This, in our opinion may be an interesting question as it is 
well known that in the presence of disorder the ballistic propagation of quasiparticles 
turns into diffusion and, in certain circumstances, into localization. It is therefore 
to be expected that entanglement itself will be affected by the presence of static 
randomness. At the level of the ground state properties the effect of disorder manifests 
in a "renormalization" of the central charge. As we will show below the dynamical 
behaviour, instead, is strikingly different in the clean and disordered cases. 
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Figure 6. The Block entropy Si for the random model for a critical value A = 1.0 
(circles) for N = 50 and m — 50. These are compared with the Refael and Moore 
prediction Eq. © (dashed line). The data have been averaged over 10 4 realizations. 



As in the clean case we analyzed the evolution of a random chain with a quench 
in the anisotropy from a non-critical A to a critical value Ai for various cases. 
Interestingly now the entropy does not grow linearly as in the non-random case. 
Although it is very difficult, in the absence of any analytic result, to ascertain the 
exact time dependence of the entropy, our data clearly indicate that the entropy grows 
logarithmically as a function of time. This is shown in Fig. where we report the t- 
DMRG results for several quenches. After a transient behaviour, all the curves in Fig. [7| 
behave like 

S t ~ nki Jt, (4) 

where k depends on the details of the quench. For example we find k ~ 0.5 for A = oo 
and Ai = 0, and k ~ 0.22 for Ao = 2 and Ai = (with In the natural logarithm). 
Clearly such logarithmic growth cannot continue indefinitely, but with t-DMRG it is 
hard to investigate times larger than those reported. 

To shed some light on the long time behaviour of Si, we consider extensively the 
diagonalization of the XX model starting from Ao = oo, as described in the appendix. 
The exact diagonalization of the XX model allows us to consider chain up to 120 sites, 
blocks of up to 50 sites and to follow the dynamical evolution to longer times as presented 
in Fig. El The analysis on the XX model was used furthermore to check the accuracy 
of the t-DMRG data. This is shown in the inset of Fig. |H] where the discrepancy is 
below 3%. 

Interestingly, Fig. |H] shows the end of the logarithmic growth of Si, even if the 
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Figure 7. Dynamics of the entropy for the random Heisenberg model obtained using 
the t-DMRG for various quenches. The solid line is a quench from A = oo to = 
for a block of £ = 10. The dotted and dashed lines are quenches from Ao = 2 to 
Ai = for two block sizes £ = 10, £ = 20 respectively. In the three cases we considered 
N = 50. Inset: the same plot but in semi-logarithmic scale. The parameters of the 
DMRG calculation are m = 60 and J St = 5 • 10 -2 . The data have been averaged over 
10 3 for the quench from A = +oo and 400 for those with A = 2. 



complete saturation is not yet reached. All the data of Fig. El are well fitted by the 
function 



An example of this fit is given in the inset of Fig. El We find that the parameter b(£) 
is slightly depending on £, being approximately 1.1 for all the cases considered. Instead 
a(£) is a decreasing functions of £, as shown in Fig. Elfor £ = 10, 20, 30, 50. We note that 
our data can be fitted by the simple scaling behaviour a(£) ~ £~ v with v = 0.16 ± 0.01, 
as shown in Fig. El At this stage however a logarithmic dependence on £ cannot be 
excluded. However, assuming that Eq. (j5J) is true for very large times, and assuming that 
the power law behaviour of a(£) is corrected, we have that for t — > oo the entanglement 
entropy saturates to u/2ln£. Thus for infinitely large times the entanglement entropy 
saturates to a value that is reminiscent of the one in the ground state, but the prefactor 
of the logarithm seems to be different (if our fit holds for large times it is exactly the 



Several considerations are in order at this stage. As compared to the clean case 
(Figs. El El EJ) the increase of the entropy as a function of time is much slower. The 
logarithmic behaviour does not follow from an extension of the argument for the clean 
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Figure 8. Dynamics of the entropy for the random Heisenberg model. Two quenches 
from Ao = +oo to Ai = for different lengths of the block: I — 30 (dashed line) and 
I = 50 (solid line). In both cases the length of the chain is N — 120 and the number 
of configurations for the averaging was 10 4 . The curves presented in this figure were 
obtained by the diagonalization of the XX model. A comparison between the t-DMRG 
and the exact diagonalization is presented in the inset for the case t = 10 and TV = 50. 



case [38] assuming that pair of particles that are emitted diffuse rather than moving 
ballistically. At this stage we cannot present arguments to derived the logarithmic 
increase of entropy in disordered chain. This behaviour is probably associated with a 
sort of entanglement localization. 

5. Conclusions 

In this paper we have presented a detailed account of the static and the dynamics of the 
block entropy for the Heisenberg model. For the static case we confirm the results of 
analytic calculations, i.e. the entropy diverges logarithmically and that the prefactor is 
proportional to the central charge c = 1 in the clean case and c e s = In 2 in the disordered 
case. For the dynamical case we have confirmed the CFT prediction that the entropy 
should grow linearly in time and that the saturation time f = !/v as a consequence of 
a reflecting wall at the boundary of the chain. The dynamical behaviour of the entropy 
after a quench depends qualitatively on the presence of disorder. Unlike for the clean 
case in the disordered case the entropy grows logarithmically with time signaling the 
possibility of entanglement localization. 
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Figure 9. Scaling of the coefficient a(£) of Eq.JSJ) with the dependence of the length 
of the block. A power law behaviour a(£) ~ l~ v (with v = 0.16 ± 0.01) fits well all 
the range considered. In the inset we show an example of a fit with the logarithmic 
behaviour suggested in Eq. (|SJ); the accuracy shown in the figure is obtained for all 
the cases analyzed in this work 
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Appendix A. 

In order to check our numerical results we compared them with the available analytic 
solution of the XX model where it reduces to a free fermion system. For completeness 
in this appendix we sketch how to calculate the entropy for the XX model with open 
boundary conditions. For A = the Hamiltonian Eq. (JI} can be rewritten, using the 
Jordan- Wigner transformation [45], as: 



N-l 




(Al) 
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where ct = (cr^ — io~ J y ). It is possible to rewrite Hamiltonian Eq. (|A.1J) in diagonal form 

w 

H = Y,E k b\b k (A.2) 
fe=i 

where the new ladder operators are connected to the old ones with an orthonormal 
transformation: 

N 
3=1 

The matrix {a k j} contains the normalized eigenvectors of the N x N adjacency 
matrix {Jkj} whose elements are defined as J\j = 2Ji<5 J) 2, J/vj = 2Jjv_i^,jv-i ^fcj = 
2J k 5 kJ+1 + 2J k _ 1 5 kJ _ 1 for k = 2, . . . , N — 1. 

To calculate the entropy of a block of size £ we use [33]: 

i 

S = ~J2 X « l0 S2 A« + (1 - Aa) log 2 (l - A a ) (A.4) 

a=l 

where A Q are the eigenvalues of the correlation matrix Cy = (c\cj), i,j = 1,£. 

As an example we consider the time dependence of the block entropy after a quench 
from A = oo to Ax = 0. The initial state is the antiferromagnetic state: 

|Vo) = c t 1 4---cj v _ 1 |0) (A.5) 

The correlations are easily evaluated using the transformation {a k j} with the result 

(c\(t) Cj (t)} = J2 a ^ a ik'e~ liEk ~ E ' k)t (H b\by |Vo) = 
k,k' 



Yl a ik a jk >a* lk a* lk ,e l{Ek (Vol c\,c y |Vo) 



k ,k' i'j' 
N-l 

= J2H a tk a jk ,a*, k a* fk ,e-^- E ^ (A.6) 

k,k' i'=l,3 

By using Eq. (|A.6|) to calculate the eigenvalues of the correlation matrix, the block 
entropy in readily evaluated from Eq. (|A.4j) . For example when the chain is homogeneous 
Ji = J we obtain: 

r^— N 

a kj = J JT—r sin k 3 ( A - 7 ) 

V 3=1 

and the corresponding energy levels are E k = 4 J cos k n where k n = (7m) / (N + 1) , n — 
1,2,. ..,N. 
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